On giant piezoresistance effects in silicon nanowires and microwires 
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The giant piezoresistance (PZR) previously reported in silicon nanowires is experimentally inves- 
tigated in a large number of surface depleted silicon nano- and micro-structures. The resistance is 
shown to vary strongly with time due to electron and hole trapping at the sample surfaces. Im- 
portantly, this time varying resistance manifests itself as an apparent giant PZR identical to that 
reported elsewhere. By modulating the applied stress in time, the true PZR of the structures is 
found to be comparable with that of bulk silicon. 

PACS numbers: 73.50.Dn, 73.50.Gr, 73.63.Nm 



As the most well studied and commercially important 
semiconductor, reports of new physical phenomena in sil- 
icon receive much attention. A recent example is giant 
piezoresistance (PZR) where the change in resistance 
of silicon nanowires due to an applied mechanical stress 
was reported to be orders of magnitude larger than that 
of bulk silicon Q. This report is highly cited [H-IH in 
part because it may represent another example of the 
effect of size on the physical properties of an otherwise 
well characterized material 

[111, El ■ Additionally, giant 
PZR is currently seen as a potential breakthrough means 
of detecting motion in nano-electromechanical systems 
5 1 where conventional detection methods lose sensitivity 
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12j. Moreover, since mechanical stress is a key 
element for performance enhancement of microelectronic 
devices [l3j], the physical mechanism behind giant PZR 
could prove to be an important ingredient in the design 
of future nano-scale transistors. As yet there is no con- 
sensus concerning the origins ofgiant PZR, although two 
models have some support. 3 13 One [3] is based on a 
surface quantization effect predicted to occur in the first 
few silicon monolayers, while the other is based on 
a stress-induced movement of the surface Fermi level in 
partially depleted structures resulting from a change in 
surface charge. The atomic length scale of the former 
seems to be in disaccord with the typical wire diameters 
reported in the literature, which are at least several tens 
of nanometers, whereas the characteristic length scale of 
the latter is the surface depletion layer width (1 nm to 10 
/j,m, depending on the doping density). It has been noted 
that the initially reported giant PZR occurred only in 
surface depleted nanowires [lj, Il5| , and subsequent claims 
of giant PZR involve depleted structures 0, la, 0]- 

Here we show that in depleted structures resistance 
changes are dominated by electron and hole trapping at 
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FIG. 1: (a) Typical layout of nanostructures, showing symbols 
used in the text and (b) SEM image of a released 2000 nm x 
2000 nm x 30 /im microwire. 



the surface. Quite unexpectedly this dielectric relaxation 
(which is independent of applied mechanical stress), re- 
sults in apparent giant PZR signatures identical to those 
initially reported in silicon nanowires [lj]. The true PZR 
can only be measured with accuracy by modulating the 
mechanical stress in time as outlined below. In all cases 
it is comparable with that of bulk silicon Q . 

A variety of unreleased and released, n-type and p-type 
microwires, nanowires and nanoribbons were fabricated 
using a top-down approach from silicon-on-insulator 
wafers of different device layer thicknesses (h) and back- 
ground doping levels (see Fig. [T]and Table HJ). The back- 
ground doping density was chosen so that the surface de- 
pletion layer width [16| Wrj > h, thereby ensuring that 
all structures are strongly depleted. All structures were 
etched using deep reactive ion etching and either ther- 
mally activated phosphorous (n-type) or boron (p-type) 
were used as dopants for the Ohmic contacts (1 x 10 20 
cm -3 ). The metal contacts and tracks were composed 
of Ti/Au for n-type contacts and Al for p-type contacts. 
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FIG. 2: (Top frame) Measurement of lDs(t) after applying a 
source-drain voltage of Vds = 0.5 V at t — across a 200 nm 
x 2000 nm x 30 (im n-type nanoribbon and alternating the 
applied stress between MPa and -13.3 MPa. Vg was held at 
V for the duration of the measurement. The inset indicates 
the sequence of measurements used in the stress modulation 
technique, where the lines are a guide to the eye. (Bottom 
frame) Relative conductance change due solely to the applied 
stress. 



For released structures, the buried oxide (BOX) was re- 
moved using an HF (50 %) etch followed by supercritical 
CO2 drying to avoid possible wire stiction. All wires and 
ribbons are aligned with the (110) crystal direction along 
which the mechanical stress is applied using a three-point 
bending method. When resistance is measured in the 
same direction as the applied stress, X, the PZR is de- 
scribed by the longitudinal piezoresistive coefficient, 



71"; = -^T 



1 AG 

X Go 



(1) 



where AG is the change in conductance from a zero-stress 
value, G - In the (110) direction, tt ; = -32 x 10 11 Pa" 1 
for bulk n-type silicon and 717 = 72 x 10 11 Pa -1 for bulk 
p-type silicon [2j. On wafer, large area strain gauges were 
also fabricated, thereby allowing the applied stress to be 
monitored in-situ. The three-point bending setup, which 
uses a piezoelectric pusher, permits rapid and repeatable 
switching between zero-stress and applied stress regimes; 
this experimental technique is different to the usual one 
where applied stress is ramped or stepped monotonically 
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FIG. 3: (a) Successive Ids — Vds measurements of a 50 nm x 
50 nm x 1 /im p-type nanowire with Vg = V. (b) Apparent 
and true PZR extracted from the Vds = 0.5 V data points in 
(a), (c) Successive Ids — Vg measurements of a 50 nm x 50 
nm x 1 (im p-type nanowire with Vds ~ 0.5 V, with an inset 
showing the apparent changes in mobility and conductance 
(calculations described in the methods section), (d) Relative 
conductance change due solely to the applied stress of various 
samples extracted using the stress modulation technique, to- 
gether with typical values for n-type and p-type bulk silicon. 



with time [j], 0, S E3] • G was measured by monitoring 
the current, Ids, through the structures under an ap- 
plied voltage, Vds, while the silicon handle was held at 
a constant voltage back-gate voltage, Vg (see Fig. [T]). 

The need to separate the time-varying and stress- 
induced resistance changes can best be illustrated with 
the experimental data shown in Fig. [5J At t = 0, 
Vds is stepped from V to 0.5 V across a 200 nm 
x 2000 nm x 30 /im nanoribbon with n-type contacts, 
while X was alternated between MPa and -13.3 MPa. 
G reduces by 27 % during the initial 1000 s of mea- 
surement, and then increases for the remaining 7000 s, 
independent of X. Using the stress modulation tech- 
nique, sequential measurements of Ids are made at zero 
stress (times t\ and £3) and with applied stress (times 
ti and £4) (see inset). The true relative change in 
the conductance AG /Go = Ids/Ids,o due to the ap- 
plied stress at £3 is found by linear interpolation to be 

[2Insfa)-lDsfo)-I D s{U)]/[lDs(t2) + lDs(U)] which 
remains at a constant -0.54 % (see bottom frame, Fig. 
I2J over time. Using Eq. [TJ 7q = -41 x 10 11 Pa" 1 , in 
excellent agreement with the bulk value for (110) ori- 
ented n-type silicon 0. If the stress had been ramped 
linearly in time, the true PZR would have been masked 
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by the non-stress-related drift of G which is 10-100 times 
larger. Indeed, the implicit assumption when using a 
linearly-ramped stress technique is that the zero- voltage 
resistance remains constant for the entire measurement. 
As is clear from Fig. [2J this is not necessarily valid for 
depleted silicon structures that, as will be seen below, are 
sensitive to surface charging. Similar dynamic changes in 
the resistance of silicon nanowires have been previously 
reported [| i, ElrSoj] . 

Figure [3k. presents the results of a measurement de- 
signed to highlight how temporal changes in zero-stress 
resistance manifest themselves as an apparent PZR. Each 
solid line represents a single Ids ~ Vds sweep for a par- 
ticular applied stress, with Vds swept from -1 V to 1 V 
in increments of 0.1 V. Two zero-stress measurements of 
Ids were taken before and after the applied-stress mea- 
surement at each applied voltage, and averaged to give 
an accurate zero-stress reference. The applied stress was 
incremented between each Ids — Vds sweep, from -13.3 
MPa up to 13.3 MPa, including an applied stress of 
MPa. The applied-stress measurements presented in this 
figure closely resemble published experimental data pro- 
duced as evidence for giant PZR (see Fig. 2b of Ref. 
ill), where the slope of each Ids — Vds curve changes in 
step with the applied stress. However, the zero-stress 
reference measurements (indicated as diamonds for the 
first TV sweep in Fig. [3^) indicate that the true PZR 
is negligible, and that the change is due solely to a non 
stress-related, time-dependent change in wire resistance 
over the duration of the measurement. Fig. [3Jd, which 
shows the apparent PZR (open triangles, calculated us- 
ing the value of current at an applied stress of MPa; the 
black line in Fig. [3k), and the true PZR (open diamonds, 
measured with respect to the zero- voltage reference). For 
clarity, the real and apparent PZRs are only depicted for 
an applied voltage Vds = 0.5 V, but similar results are 
obtained at all voltages. The apparent PZR has an expo- 
nential dependence on X like the curves obtained in Ref. 
[l| (see Fig. 2c of that article) although other forms (in- 
cluding those labelled C, Z and I in Ref. 1) are observed 
at different moments along relaxation curves of the type 
shown in Fig. [2] Similarly, changes that are even larger 
or of opposite sign are equally possible. 717 obtained from 
the apparent PZR in Fig. [3b is 717 = 450 x 10 11 Pa" 1 
around X = MPa, and without use of the stress modu- 
lation technique, this value is indistinguishable from the 
true PZR. 

By measuring the change in Ids in response to a 
change in Vg, and by assuming a linear relationship be- 
tween the slope of the Ids — Vg characteristic and the 
mobility, the giant PZR was attributed to a mobility 
variation. [l[ . This measurement is replicated using a 50 
nm x 50 nm x 1 /im p-type nanowire. The results are 
presented in Fig. [3fc, demonstrating that the apparent 
mobility extracted from the slope of an Ids ~ Vg mea- 
surement can also change over time independently of the 
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FIG. 4: Measurements and modelling of the change in AVa 
due to oxide charge trapping and the values of Ios for a 200 
nm x 2000 nm x 30 /im n-type nanoribbon, measured in 
% and 40 % relative humidity. The source-drain voltage Vds 
was held constant at 0.2 V, while step changes were applied 
to the gate voltage Vg as indicated in the top frame. 



applied stress. The sign and magnitude of this change is 
the same as the apparent stress-dependent conductance 
change. The true stress-dependent conductance change 
is consistent with bulk silicon PZR (see inset). Figure [3Ji 
shows the true values of AG/Gq as a function of stress 
obtained for four different, depleted silicon structures, 
together with the values expected from bulk silicon. Re- 
gardless of lateral wire size, or whether the device is re- 
leased or not, the true PZR compares well with that of 
non-depleted, bulk silicon. Results from all measured 
samples can be found in Table HI 

To better understand the origin of the dynamic con- 
ductance changes, a similar approach to that used in a 
study of oxide traps in MOSFETs is employed, in which 
a measurement of the subthreshold current, 



Ids = hexp[{V G + AV G )/S], 



(2) 



is made.[2l[ Here Iq is a constant, S is the sub-threshold 
slope, and AVg is the shift in the effective gate voltage 
due to trapped charge. AVg a t time t can then be ex- 
pressed as 



AV G {t) - AV" G (0) = Sin 



Ipsit) 
Ids(0) 



V G (t) + V G (0), (3) 
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where S can be measured directly from a rapid Ids ver- 
sus Vc measurement. Measured Vait) — Vg(0) data can 
then be fitted using a semi-empirical model for positive 
and negative charge trapping in an oxide which is based 
on an existing model [22| for positive charge trapping 
in SiC>2 layers of metal-oxide-semiconductor (MOS) ca- 
pacitors. The model assumes that the sheet density N 
of trapped charge has a saturation value that de- 
pends linearly on the electric field across the oxide £ ox , 
so that Nrx, — a^ ox where a is a constant with units of 
cm _1 V w 1 . N approaches Noo with a time constant r, 
so that dN/dt = -(N(t) - N^/t, where r is related 
to the capture cross-section of the trap and the current 
density through the oxide [H, [23| . The contribution of 
the trapped charge to the gate voltage is given by [22[ 
AVq = qNd ox /e ox , where d ox is the oxide thickness, q is 
the charge on an electron, and e ox is the oxide permittiv- 
ity. Figure [4] presents measured values of Vg (t) — Vq (0) 
for a 200 nm x 2000 nm x 30 /im n-type nanoribbon, 
where Vds = 0.2 V and a series of positive and negative 
steps were applied to Vg ■ Two sets of measurements were 
taken: one in a dry nitrogen atmosphere (0 % relative hu- 
midity) and one in a relative humidity of 40%. Each set 
of data is modelled using one type of electron (hole) trap 
characterised by N e (0), a e and r e (Af^(0), ah and t^). 
Excellent agreement is reached between the model and 
the measured values of Vc(t) — Vg{0) over each 20- hour 
measurement period. The calculated current is deter- 
mined by using the modelled values of Vc(i) in Eq. [2] and 
by adjusting the parameter J . This closely matches the 
measured current as shown in the bottom frame of Fig. |4] 
This treatment yields several important insights into the 
origin of the dynamic conductance changes: i) The ob- 
served changes in current are the result of charge trapping 
of electrons and holes in SiC>2 layers at the wires surfaces; 
ii) The changes in current are much more rapid at higher 
relative humidity, consistent with the presence of water- 
related charge traps [22|, |24[ ; iii) The fitted values of the 
constant ah for the hole traps (1.2 x 10 12 cm _1 MV _1 at 
% relative humidity and 1.5 x lO^cm^MV -1 at 40 
% relative humidity) are almost identical to values ob- 
tained in MOS capacitor oxides 22]; iv) Electron traps 
are associated with the rapid initial change in Ids, while 
hole traps are associated with the slower change in the 
opposite direction, consistent with observations of charge 
trapping in MOS capacitors [23J. This is strong evidence 
that the observed dynamic conductance changes are due 
to water-related charge traps in the oxide layer at the 
silicon surface. Consistent with this, the apparent giant 
PZR of Ref. [H was also shown to depend strongly on the 
characteristics of this oxide layer. 

In conclusion, charge trapping and de-trapping can 
mask the true PZR of depleted structures and mistak- 
enly lead to claims of an apparent giant PZR. In more 
than 20 different surface depleted samples, the true PZR 
is found to be comparable with known values for bulk sili- 



TABLE I: True 7r; measured on all samples. 
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con (see Table U|).[25j This is in stark contrast to previous 
reports of giant PZR in structures of similar dimensions 
and doping levels. While this does not rule out giant 
PZR in depleted silicon structures, future claims must 
conclusively demonstrate that any measured resistance 
change be solely due to the applied stress. 
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